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Abstract. Let be the discrete Laplace operator acting on functions (or rational matrices) 
/ : Ql Q, where Qz, is the two dimensional lattice of size L embedded in 1,2 ■ Consider a 
rational L x L matrix H, whose inner entries Tiij satisfy V^Tiij = 0. The matrix H is thus the 
classical finite difference five-points approximation of the Laplace operator in two variables. 
We give a constructive proof that Ti is the restriction to Q_t of a discrete harmonic polynomial 
in two variables for any L > 2. This result proves a conjecture formulated in the context of 
deterministic fixed-energy sandpile models in statistical mechanics. 
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Introduction and Motivation 

An interesting class Ml of L x L matrices Tl with rational entries and a 
related vector space of polynomials in two variables arise in some theoretical 
physics models, the so-called deterministic fixed-energy sandpiles (DFES) with 
Bak-Tang-Wiesenfeld (BTW) toppling rule [3j. 

Introduced for the first time in [4J by imposing a global energy conservation 
constraint on its dissipative counterpart [2j, DFES is a deterministic cellular au- 
tomaton, in which two-dimensional configurations (represented by square ma- 
trices with integer elements Zij{t)) evolve in discrete time steps t according to 
a precise parallel updating rule. 

The main feature of DFES is that, in contrast with the dissipative model, 
only a small part of an a priori huge configuration space is dynamically explored, 
and the system enters a periodic orbit after a surprisingly short transient. This 
is a clear indication of the existence of many hidden conservation laws (HCL) 
which split the wide configuration space into dynamically intransitive, and thus 
much smaller subspaces [3]. 

Few of those HCL were identified in a non-systematic way in [1] and can be 
represented in the form: 



^L[/](t) 



mod L 



where the sum runs over the integer coordinates of the two-dimensional L x 
L lattice sites, Zij{t) is the integer value taken by the entry at time t 

and f{i,j) is a L x L matrix with rational entries. The interest is then in 
characterizing the generating functions (GF) of HCL, i.e. the class of inequivalent 
matrices / such that $L[/](t) is a HCL (^>L[/](t + 1) = '^L[f]it) for all t). 

Bagnoli et al. [1] gave the following three GF: fi = i, f2 = j and /s = 
i'^ — j^. An intriguing observation is that, when thought as functions on the 
whole {f{x, y) :M? M), those three GF belong to a special vector space of 
polynomials in two variables, which we call discrete harmonic polynomials (see 
Def. [5]). It is then appealing to conjecture that this should be a general feature 
of any GF of a HCL. 

In fact, an exhaustive characterization of GF has been given in [3j from a 
completely different perspective, i.e. without any reference to polynomials, but 
working simply on the matrix representation of those GF. 

It was proven in [3] that a functional of the form ([T) is a HCL if and only if 
its GF is a inner-harmonic matrix of size L (see Def. 



^In appendix B of [3] the necessary and sufficient condition is expressed in terms of K- 
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The purpose of this paper is to provide a rigorous link between the exact 
characterization of HCL in terms of matrices [3] and the conjectured polynomial 
form for any GF. More precisely, we will prove that every inner-harmonic matrix 
of size L (i.e. any GF of a HCL in the sandpile context) can be represented (non 
uniquely) as the restriction to the two-dimensional discrete lattice of a discrete 
harmonic polynomial in two variables . 

The paper is organized as follows. In Section [T] we set up notations and ba- 
sic definitions, providing in particular the notions of i) inner-harmonic matrix 
of size L (Dei. |4|), in terms of the well-known five-points formula for the dis- 
cretization of the Laplace operator on a 2d lattice, and ii) discrete harmonic 
polynomial (Def . [5]) . In Section [21 we enunciate the main theorem and provide 
the algorithmic procedure for finding the discrete harmonic polynomial which 
interpolates any given inner-harmonic matrix of size L > 3. In the same section, 
we provide a stepwise example of application, together with pointers to subse- 
quent lemmas needed for the proof. Section [3] is devoted to conclusive remarks 
and hints for future works, while a basis of discrete harmonic polynomials up 
to degree 9 is given in the Appendix. 

1 Definitions 

We define Ql as the two dimensional lattice embedded in Z2, i.e.: 

QL = {{i,j) eZ2\0<i,j <L-1} (2) 

Definition 1. The inner sublattice of is the set: 

Q[ = {{^,j) (^QL\l<i,j <L-2} (3) 

The discrete Laplace operator is defined as the classical finite difference 
five-points second order formula for the approximation of the Laplace operator: 

Definition 2. Let / : Ql Q- The discrete laplace operator acts on / 

as: 

= ^f{iJ) - fii - 1,J) - f{i + - /(^,J - 1) - f{iJ + 1) (4) 
where p := /l^t^. 

The generalization to functions / : ^ R is straightforward (consider 
/ = /Mn this case). 

harmonicity, and strictly speaking this is not equivalent to inner- harmonicity. However, it 
can be proved that for every K-harmonic function there exists an inner-harmonic function 
which belongs to the same equivalence class, i.e. generates an equivalent HCL. Thus, it is not 
restrictive to work with inner-harmonic matrices, as we will do from now on. 
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Definition 3. Let Ml be the set of rational Lx L matrices and ^ = {f\f : 
Ql Q}- We define the invertible map ■ ^ ^ Ml {L- correspondence) 
through the following: 

^(/i) := n (5) 

where h{—l + j, —i + L) := l~Li,j- 

Through \1', the lower left corner of V. is mapped to the point (0, 0). 

Definition 4. A L x L rational matrix Hi is called inner-harmonic matrix 
of size L (L > 2) if the following property holds {hi = '^~^{Hi)): 

(vVi)(i,j) = (6) 
as in the following example, where we restrict for simplicity to integer entries: 



f 2 








1 





1 


2 \ 





2 


1 


2 





2 


1 


1 


7 





6 


-4 


6 


2 


1 


25 


-14 


26 


-28 


24 


1 


2 


106 


-107 


140 


-158 


117 


2 


2 


504 


-660 


799 


-861 


600 


1 


V 1 


2568 


-3836 


4577 


-4685 


3143 





Definition 5. A polynomial P{x, y) is called discrete harmonic polynomial 
if (V2p)(a;,y) = ^{x,y)e^^. 

Examples of discrete harmonic polynomials are Pi{x, y) = x^—y^, P2{x, y) = 
x^ - 3xy'^,P3{x,y) = xy. 

The set of discrete harmonic polynomials of degree g will be denoted as B*. 

Definition 6. We say that a polynomial P{x, y) interpolates a LxL matrix 
H if P{i,j) = h{i,j), where *(^) = H. In this case, we write P = H. 
Note that: 

Remark 1. Discrete harmonic polynomials are generally not harmonic in 
M^, i.e. solutions of the continuum Laplace equation V^P = 0. Generally speak- 
ing, every polynomial V{x,y) in two variables belongs to one of the following 
classes: 

• V{x,y) is neither harmonic nor discrete harmonic. Example: V{x,y) = 
x^ + y^ 

• V{x, y) is harmonic but not discrete harmonic. Example: 'P{x, y) = x'^ — 
dx'^y'^ + 

• V{x, y) is discrete harmonic but not harmonic. Example: 'P{x, y) = x'^ — 
2x^ — 6x^y^ + y^ 
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'P{x, y) is both harmonic and discrete harmonic. Example: 'P{x, y) = xy 



Remark 2. Given a discrete harmonic polynomial P{x,y), it obviously in- 
terpolates an inner-harmonic matrix TCl on Qj,. For example, the polynomial 



P{x,y) 



X 



3xy interpolates the following matrix on Q 



7- 



/ 216 


198 


144 


54 


-72 


-234 


-432 \ 


125 


110 


65 


-10 


-115 


-250 


-415 


64 


52 


16 


-44 


-128 


-236 


-368 


27 


18 


-9 


-54 


-117 


-198 


-297 


8 


2 


-16 


-46 


-88 


-142 


-208 


1 


-2 


-11 


-26 


-47 


-74 


-107 


^ 

















) 



(8) 



The converse is not trivial for any L > 2: while it is straightforward to find an 
interpolating polynomial ^{x,y) for any given inner-harmonic matrix through 
any of the known Polynomial Interpolation formulas in two variables [5j, the 
resulting $ is generally not discrete harmonic in M? (and incidentally neither 
harmonic). This can be seen easily by referring to the widely used Bilinear 
Interpolation formula (see e.g. [8]), the extension to the two-dimensional lattice 
of the well-known Lagrange interpolation formula in Id: 



L-l 



^{x,y) = J2^hkYl 



h,k 



j=0 



. L-l 

-TT- 

r— 



(9) 



where Zhk = "Hh^ki the sum runs over the sites of the matrix and the products 
over rows and columns respectively. Note that deg(<I>) = 2(L — 1). 

It is then possible to interpolate the following simple inner-harmonic matrix 
of size L = 4: 

/ 27 18 -9 -54 \ 
8 2 -16 -46 
1 -2 -11 -26 
-3 



7^4 



V 



(10) 



The bilinear interpolating polynomial is the following: 

11 o 1 , 11 121 5 
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^hA^i 2/) = -3 + y X - 3x^ + -x^ + —y - -^xy + -x^y - -^x^y+ 

- 3?/^ + -^xy^ 



3x^2/2 + ia;3y2 ^ 3^3 



11 Q loo 

— xy'^ + -x^y^ 

12 2 



x^y^ 
(ID 



and a straightforward calculation yields V'^{^f{^){x,y) 7^ in 
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In the following section, we shall provide the enunciation of the main result, 
a stepwise example of application of the algorithm, and a constructive proof of 
the main theorem. 



2 Interpolation by discrete harmonic polynomials: 
main result and algorithm 

We enunciate our main result: 

Theorem 1. Let 7i be an inner-harmonic matrix of size L > 2. There exists 
a discrete harmonic polynomial P{x, y) of degree less or equal to 2{L — 1) such 
that P interpolates TL on Q/,. 

Before getting to the technical points, it is informative to provide an example 
of how our algorithmic procedure roughly works. 

Let us consider the inner-harmonic matrix Ti := T~Li in (jlOp . 



First step: 



3) minor 








8 2 


-16 




1 -2 


-11 




-3 






(12) 



Second step: 

Apply Lemma [3] and find a discrete harmonic polynomial! P'^^\x,y) = TL^^^: 

, N 15 1 o 3 Q 1/1 15 

P^i) (x, y) = -3 + — a; - -x^ - -x^ + -x^ + — y+ 



Third step: 

Evaluate P^-^\x,y) on the lattice Ql = Q4, obtaining the matrix 7Y4: 



7^4 



/ 24 |l8 

8 2 

1 -2 

\ -3 



-9 
-16 
-11 



-57 \ 
-46 
|26] 
-3 / 



(14) 



^Note that this polynomial does NOT coincide with the bilinear interpolating polynomial 
we would get for the same matrix. 
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Note that i) 7^4 / TL ii) the sites (1,3) = 18 and (3,1) = —26 are uniquely 
determined by the discrete harmonicity requirement and thus coincide in the 
two matrices. 

Fourth step: 

In order to amend the other mismatching entries along the border, compute the 
four (L)-Polynomials (Lemma H]): 



ei(x,y) = 384x 



65x'' 



3x^ + 



5167/ + 332xy+ 



+ 465x^7/ - Amx-^y + 105a;*y - Qx^y + 776y^ - 1095x?/"'+ 



+ 360x^7/2 + 30x^7/2 



15x'^y2 



225y^ + 460x7/2 



210x^7/3+ 



+ 20x^7/2 



55y^ 



15x/ + 15xS'' + 21y^ 



6(x, y) = 240x - 386x2 _^ -^35^3 ^ 25x'^ 



15x^ + x^ 



6X7/ — 7/ 

- 168y - 152xy+ 



(15) 



+ 555x2y - 280x^7/ + 15x''y + 6x^7/ + 3267/^ - 255xy2+ 

1957/2 + 260xy2 
75x7/"' + 15x^7/^ + 37/^ + 6X7/ 



180x^^2 + 150x^7/2 



30x2y3 + 



20x^y^ + 357/^ 



(16) 



^3(x, y) = 516x - 776x2 _^ 225x3 + 55x'' - 21x^ + x*^ - 3487/ - 332x7/+ 



+ 1095x2y - 460x^7/ + 15x'^y + 6x^y + 656y2 



465x7/"^+ 



360x27/2 + 2 10x^7/2 



15x^y2 



375y3 + 440x7/3 



30x2y3+ 



20x^7/3 + 65y^ 



105x/ + 15x2/ + 3y^ + 6xy^ 



(17) 



^4(x, 7/) = 1644x - 2852x2 + 1305x3 - 35x'^ - 69x^ + 7x^+ 
- 1644y + 3225x2y - 2130x3y + 345x*^y + 2852y2 + 



3225xy2 



690x3y2 



105x^y2 



1305y3 + 2130xy3+ 



690x2y3 + 35/ 



345x/ + 105x2y4 + 69y^ 



7y^ 



(18) 



Those have the remarkable properties to be i) discrete harmonic in M2 ii) 
almost everywhere on Q^, except one single entry (two for ^4). In particular, 
^k=ik, where: 

/ 71 \ 


V / 



6 



(19) 



/ 72 \ 





V / 



(20) 
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/ \ 





V 73 / 

/ 74 \ 
-74 


V y 

where, for the particular choice of the basis polynomials used to build up the 
^{x,y) (see LemmalUfor details), we have (71,72,73,74) = (—720, —720, 720, —720). 



^4 



(21) 



(22) 



Fifth step: 

Define the sought interpolating polynomial P{x, y) for 7i as: 



P{x,y) = P^^Hx,y) + ^ZkUx,y) 



k=l 



where Zk are parameters to be determined, and compute P{x,y) on Q/,: 



V 



/ 24 - 720zi 18 -9 - 720^4 -57 - 720^2 \ 

8 2 -16 -46 + 720^4 

1 -2 -11 -26 

\ -3 -3 + 72OZ3 J 



(23) 



(24) 



Sixth step: 

Compute {zi, Z2,Z3, Z4) by requiring V = 7i: 



24 - 720zi 
-57 - 720z2 
-3 + 720^3 
-9 - 720^4 



which gives: 



27 



-54 



zi = -1/240 

Z2 = -1/240 

Z3 =1/240 

Z4 =0 



(25) 
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Substituting (p5|) back into (|23|) . the final result is obtained: 

, ^ 69 59 2 31 o 25 4 1 . 1 e 103 

J-^(x,y] = —6 -\ X -\ X X -\ X X X -\ y+ 

^ 20 60 16 48 80 240 20 

533 To 13 s 7 4 1 . 119 2 

xy X y H x y x y H x y y + 

60 ^ 16 ^ 12 ^ 16 ^ 40 ^ 60 ^ 

+ TT:xy^ - 3x^y^ + -x''y' + —x^y' + —y' - -xy" + -x^y'+ 
16 8 16 16 6 8 

+ — x^y^ + — / - —xy^ - —x'^y'^ - —y^ + —xy^ + — / (26) 
12 ^ 48 16 16 80 40 240 ^ ' 



Note the difference between (j26p and (jlip although they interpolate the very 
same matrix (llOp . The degree of P is 6 = 2(L — 1) as stated in Theorem [TJ 

This procedure can be iterated without difficulties up to interpolating inner- 
harmonic matrices of any size through a repeated application of Lemma [5l 

Hereafter we shall provide several preliminary lemmas which are essential 
for the proof of the main result and have been hinted previously. 

Lemma 1. Let k > 0. Then is a vector space of dimension 2. 

First, we easily prove the following statement: let P^v be the set of two 
variables polynomials up to degree N and let Pi{x,y) € Pat. Then P2{x,y) = 

V^Pi{x,y) eFN-2- 

In fact, we notice that the following properties hold: 

V2(ax" + by"') = aV^x-^) + bV^{y"'') Linearity (27) 

V2(a;"y'") = x^V^iy"") + y'^V2(x") Leibniz rule (28) 

Furthermore, for every one- variable monomial in x (or y), it is straightfor- 
ward to prove the following: 

V2,» - / -2 Elr„="= .f n - even 

"l-2Et-„'"n.i,)-'=« .f—dd '''' 

Therefore, applying the Laplace operator to a one-variable monomial of de- 
gree n, we obtain a linear combination of one-variables monomials up to degree 
n — 2. Thanks to ()27p and ()28p . we can conclude that the same holds also for 
two- variables polynomials. [QED] 

It is well-known that Pat is a linear vector space, with dim(P7v) = J2n=o(.^~^ 
-^■j _ iN+i)(N+2) ^ According to the previous results, we call Hn '■ Pa^ — > ^N~2 
the following linear map: 

UMiP{x,y)) = {S/^P)ix,y) (30) 
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Then, we call Bat = ker(n7v), i.e. the following vector subspace of Pat: 

Djv = {P(a;,y) GPjv: (V2p)(x,y) = 'i{x,y) eR''} (31) 

The elements of Djv are discrete harmonic polynomials. The dimension of Bjv 
can be found simply applying the Rank- nullity theorem to the map XIat: 



dim(Djv) = dim(Pjv) - dim(Pjv_2) 



{N + l){N + 2) N{N-1) 



Let be the following vector subspace of B^r: 

Bjfc = {P{x, y) G BatI P's degree is exactly k < N} 



2N + 1 
(32) 

(33) 



Obviously, dim(D*) = dim(Dfc) - dim(Dfc_i) = 2. 

Therefore, for A; > we can always find two (and not more) linearly indepen- 
dent discrete harmonic polynomials, i.e. elements of Bat, with the same degree 
k. 

Following the standard algebraic procedure, it is quite easy to build up a 
complete basis B* = {e*, . . . jCgAr+i} for Bat, starting from the canonical basis 
in Pat: 

Bn = {l,x,y,x'^,xy,y'^,...,y^} 

Throughout this paper, we will refer to the basis {Uk{x,y)} listed in the 
Appendix. 

Lemma 2. For every square matrix with an arbitrary fixed rational contour, 
there exists one and only one inner-harmonic completion. 

Let FQ{i,j) : Ql Q and let its (4L — 4) border sites be forced to assume 
rational values belonging to the set fl. 
In matrix form, we have: 



( Zi Z2 

Z4L-4 Xl X2 

Z4L-5 XL-l XL 

\ ZzL-2 Z3L-3 ■ ■ ■ 



ZL \ 
XL-2 ZL+l 
X2(L-2) ZL+2 

Z2L Z2L-1 J 



(34) 



The nested (L — 2) x (L — 2) submatrix has unknown entries xj € Q. 

We prove that, for each set O, there exists one and only one submatrix F^ 
with rational entries such that Fn(i,j) is inner-harmonic. 
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If we impose the inner-harmonicity condition on Fq, we get the linear system 
Ax = f]{{zk}), where A is the following (L — 2)^ x (L — 2)^ matrix: 



41 - H = 41 





I 


••• 


\ 


I 


H 


I 


... 





I 


H I 


... 








I 


V 






I h/ 



(35) 



I is the identity matrix (L— 2)^ x (L— 2)^, I is the identity matrix {L — 2) x (L— 2), 
is the null matrix (L — 2) x (L — 2) and H is a well-known matrix describing 
the Hamiltonian of nearest-neighbor hopping on a one-dimensional lattice (see 
[7] and references therein): 



H 



/ 1 
10 10 
10 1 



\ 



\ 





1 



(36) 



The vector ff depends on the fixed contour values. Its entries are of the 
following forms: 



Za + Zj3 if Xj is a corner site of 







if Xj is a border site of F^, but not a corner site 
otherwise 



Since the matrix (I35p is diagonal predominant [BJ, the system admits one 
and only one solution in Q. [QED] 

Corollary 1. If Q = {0, . . . , 0}, then F^ is the null matrix L x L. 

Lemma 3. Given a 3 x 3 inner-harmonic matrix A, it is always possible to 
find a discrete harmonic polynomial P{x, y) with rational coefficients and degree 
4 such that P = A on Q3. 

For L = 3, there are 8 sites along the contour. Choose the following set of 
discrete harmonic polynomial^ (see Appendix): 



{Uo{x,y),...,U6{x,y)}U{U8{x,y)} 



(37) 



Obviously, infinitely many other choices are equally possible. 
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The sought polynomial P{x, y) satisfying the Lemma may be written as a linear 
combination of the polynomials in (I37p . with unknown coefficients aj [j = 
1,...,8). 

The condition that P = A translates into a linear system with 8 equations 
in the unknowns Oj, whose matrix of coefficient for the choice ()37p is: 



M 



/ 1 




















\ 


1 





1 





1 





1 


-1 


1 





2 





4 





8 


8 


1 


1 








-1 


1 





1 


1 


1 


2 


2 


3 


-11 


2 


-15 


1 


2 








-4 


8 





16 


1 


2 


1 


2 


-3 


2 


-11 


-9 


V 1 


2 


2 


4 





-16 


-16 


-72 / 



(38) 



The determinant of M is nonzero. Thus the polynomial interpolating the contour 
(and for Lemma [2] also the central site) always exists and has degree 4. [QED] 

Lemma 4. For every L > 3, there exist four discrete harmonic (L)-polynomials 
^2, £.3 and ^4, whose degree is less or equal to 2L, such that = Z'^^\ k = 
1,2,3,4. The entries of the matrices Z^^\ k = 1,2,3,4 are all except: 

(1) the entry (0, L) for Z^^); 

(2) the entry {L,L) for Z^^); 

(3) the entry (L,0) for Z^^),- 

(4) the entries (L — 1,L) and {L,L — 1) for Z^^^ ; 

As it was evident from the example of application, the (L)-polynomials have 
the following effect. Given a (L + 1) x (L + 1) inner-harmonic matrix Q and a 
discrete harmonic polynomial P{x, y) interpolating the lower-left minor (L x L) 
of Q, those polynomials neutralize the mismatch between 4 sites along the border 
of G and the values assumed by P{x,y) on Ql+i- 

We prove now the existence of ^i{x,y). For the others, the procedure is 
completely analogue. 

Consider a set of 4L linearly independent discrete harmonic polynomials 
Pk,six,y), where k = 1, ...,2L is the degree, and s = 1,2 which do not contain 
the constant term. 

We write the sought ^i(x, y) in the form 



Ciix,y) = aiPi^i{x,y) + a2Pi,2{x,y) + ... + a4LP2L,2{x,y) (39) 
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To determine the 4L unknowns aj, we require that ^i{x,y) should be zero 
on i) the border sites of the lower-left minor A4 of Z^^^ ii) four other points in 
Z2, precisely: (L - 1, L), (L, L), (L, 0), (L + 1, L). 

This translates into a linear homogeneous system S in 4L equations for the 
4L unknowns ai, 04/,. Note that the site (L, L — 1) is automatically zero due 
to the harmonicity condition. 

The first row of the matrix of coefficients for S is given by: 

^1,1(0, 0), Pl,2(0, 0), P2L,2(0, 0) (40) 

and these values are all zero because the polynomials do not contain the constant 
term. 

Thus, the determinant is zero and the homogeneous system has an infinite 
non-zero solutions set {ai, ...,04^}. Since the polynomials Pk,six,y) are linearly 
independent, the obtained polynomial cannot be identically zero by definition. 

Now, let ^1 be defined by a non-zero solution {ai, ...,0^1}- Being zero along 
the contour of , it is zero inside M because of the Corollary [TJ 

It is also zero by the discrete harmonicity relation on sites {k,L), k = 
1, .., L + 1, and sites (L, k), k = 0, L + 1. Instead, it is required to be nonzero 
on the site (0, L). Indeed, we can prove that this is the case by contradiction. 
Assume that Ci(0, L) = 0. We have: 

Ci{j,k)=0, k = L, i = 0,1,...,L + 1 . (41) 

Let: 

n = L/2 + 1 m = L/2 for even L (42) 

n = {L + l)/2 + l m = {L + l)/2-l for odd L (43) 

Due to the harmonicity relation, there is an integer J, 0< J < L + 1, such 
that: 

^i( J, L — 1 -|- i) = for every i = 1, n 

^i(J, — A;)=0 for every A; = 1, m 

This means that the one variable polynomial ri{y) = ^,i{J,y) has 2L + 1 
zeros: but this is absurd, since its degree in y is at most 2L. Therefore ^i(0, L) 7^ 
0.[QED] 

Lemma 5. Let A be a inner-harmonic matrix of order L, and A' the {L — 
1) X (L — 1) lower-left inner-harmonic minor of A. Let xi^j v) 0. discrete 
harmonic polynomial of degree h interpolating A! . Then, it is possible to de- 
fine a discrete harmonic polynomial a{x,y), of degree k = max[2(L — 1), h], 
interpolating A. 
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Define: 



f 

Si 


= Site (0, L 


-1) 




S2 


= Site (L - 


2,L- 


1) 


S3 


= Site (L - 


1,L- 


1) 


S4 


= Site (L - 


1,L- 


2) 


.55 


= Site (L - 


1,0) 





and denote xisk) '■= Xk and A{sk) ■= flfc for simplicity. 
We write the sought (t{x, y) in the form: 



o-(.x, y) = x{x, y) + X] ^kik{x, y) 



(44) 



k=l 



where the Ckix,y) are (L-l)-polynomials as defined in Lemma U and are 
coefficients to be determined. The degree of each of the (^k is at most 2(L — 1), 
confirming the statement of the Lemma about the degree of a. 

We note that a{x,y) = x{x,y) on the sites of A', since all the (L-1)- 
polynomials assume value there. 

The values assumed by the polynomial x on the North and East borders of 
Ql are uniquely constrained by the harmonicity condition, except the five sites 
Sk- In general, at ^ Xk- 

For example, for L = 5 we have the following schematic situation (compare 
with ([Ml)): 

/ □ ■ ■ □ \ 



A 







□ / 



(45) 



where: 



□ 





Sites in A' , where x = A 

Sites where x = ^ by harmonicity 

Sites (51,53,55) where x^^A 

Sites (52,54) where x^^A, but mutually constrained by harmonicity 



Indeed, the discrete harmonicity condition, applied to A and X) requires that: 



X2 + X4 = 02 + 04 



(46) 
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Given that Ckisi) = for /c 7^ 1 and ^i(si) = 7 7^ (Lemma H]), we get from 
equation (fill) : 

'^(•si) = Xi + 2:17 = ai (47) 

This determines zi as zi = (ai — Xi)/7- 

The same procedure applies to the sites S3 and S5, determining Z2 and z^: 
note the shift of indices, reflecting the fact that we have five sites and only four 
(L-l)-polynomials. 

In fact, the polynomial ^4 has to be nonzero simultaneously on both sites 
S2 and S4, and by harmonicity ^4(52) = —(,4{si). This constraint, however, is 
compatible with the correct definition of 24 and therefore of a{x,y). 

Indeed, define uj = ^4(52) = —^4(54). Equation (l44l) requires evidently that 
o'(s2) = X2 + z^oj = 02 and £1(54) = Xi — z^to = 04. Both equations are obviously 
satisfied by Z4 = (02 — X2)/^ thanks to ([Ml). 

Thus, the coefficients zi, .., 2:4 in (fHI) are uniquely determined and the poly- 
nomial a(x, y) interpolating A exists. [QED] 

We are now able to provide a proof of Theorem [TJ 

We only need an iterative (or "telescopic") application of previous results: 
starting from TY, we drop the upper row and last column on the right, defining 
the minor 

If we can find a discrete harmonic polynomial xi^iV) = 'H^^\ such that 
deg(x) < 2(L — 1), the Theorem follows via Lemma [5l otherwise, we drop the 
upper row and last column on the right of TL^^^ again, and restart the procedure. 

This process is consistent, because the minors iteratively defined continue 
to be inner-harmonic. 

Suppose that we have finally found the minor (whose size is L — n) of 
H^"'~^\ admitting an interpolating polynomial X(n)(3^iy) such that deg(x(n)) < 
2(L — 1). By Lemma[5l the minor 7i^'^~^^ can be interpolated, and so on, up to 
interpolating 7i. 

Since at least for L = 3 the interpolating polynomial always exists (Lemma 
[3|), in the worst possible case the telescopic algorithm will start from the 3x3 
lower left minor of Ti. , and will eventually produce the desired result by repeated 
applications of Lemma [5l [QED] 

3 Final remarks 

In this note, we have developed a "telescopic" technique to interpolate an 
inner-harmonic matrix of size L by a discrete harmonic polynomial of degree 
less or equal to 2(L — 1). 
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The solution we have presented proves a conjecture about hidden conser- 
vation laws in the context of some statistical mechanics models, namely the so 
called fixed-energy sandpiles with deterministic BTW toppling rule. 

We remark that the algorithmic procedure we devised should be regarded as 
a mere tool for the proof, and by no means is meant to provide a computationally 
efficient and robust interpolator for inner-harmonic matrices. 

As a final point, we wish to give here a short survey on other related questions 
and problems which have not been addressed in this paper and could be worthy 
of further investigations. 

(1) Discrete harmonic polynomials of minimal degree: the constructive proce- 
dure outlined in section [2] does not lead to an uniquely defined interpolat- 
ing polynomial. A natural question to ask is what the minimal attainable 
degree of such a polynomial is, and how to build it up. 

(2) A related combinatorial problem: Another class of matrices {M^) with 
integer entries and closely related to A^^, emerges in [3] and proves to be 
connected to deep symmetries of the evolving rule of that model. 

The main features of Al^ are: 

• Condition ([6j) holds modulus the size L of the matrix. 

• Cyclical border conditions are imposed and condition Q holds for 
border sites as well. 

• Entries are bounded by an integer M. 

An interesting problem in analytical combinatorics, with many possible 
consequences on the underlying physical issue, is to count the number of 
those matrices for fixed L and M. 
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Appendix: list of Discrete Harmonic Polynomials 

We report here a basis of discrete harmonic polynomials up to degree 9 that 
we used repeatedly throughout the paper: 



Uo 
Ui 

U2 
U3 
Ua 
U5 
Ue 
U7 
Us 
U9 
Uio 

Uu 
U12 

Uu 

Ui5 

C/16 

Ui7 



u 
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x,y 
x,y 

x,y 
x,y 
x,y 
x,y 
x,y 
x,y 
x,y 
x,y 
x,y 

x,y 

x,y 
x,y, 
x,y 

x,y 
x,y 

x,y 

x,y 



= 1 



X 

xy 

2 2 

X -y 



-3x^y + y^ 
3xy^ 



X 



x^y — xy^ 
x^ - 2x^ - 



6x y +y 



5x^y - lOx^y-^ - lOx^y + y'' 



X 



x^y 



10 



10 



Wxy^ 
-xy^ + xy^ 



= -Ibx'^y^ - lOx"^ + lOx^ + 15x\^ + SOx'^y^ - y^ + 

= 35xS^ + lOx'^y - 21x'^y^ - 70x'^y^ - lOx'^y + - 7x^y 

= -2lx^y^ - TOx^y^ + 35x^/ - 7xy^ + 70a;/ - 70xy'^ + 

= -7x^y^ + Tx'^y^ - y-'^^'^y'^ - y + ^v'^ + 14x?/^ + x'^y 

= -140x^y2 + TOx"^/ - UOx'^ + 166x^ - 28x'^y^ + 280x^y'^+ 
+ 560x^^2 + y8 _ 28/ + x^ - 28x^y^ 

= 126a;^/ - 252x^y^ - 84x^y^ - SAOx^y^ + 840x^/ + 9x/+ 
- 252x/ + 1260.xy^ - 1026xy^ + x^ - 36x'^y^ 



= SAOx^y^ + UQx^y^ + 1260x^y - 252x2/ - SGx^'^ - 840x2y^+ 
- 1026x2 y + / + 9a;^y - 84x^/ - 252x^y 



(48) 
(49) 
(50) 

(51) 
(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 

(59) 

(60) 
(61) 
(62) 

(63) 

(64) 

(65) 

+ 

(66) 
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